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Abstract 

The cortical dynamics of somatosensory processing can be investigated using vibrotactile 
psychophysics. It has been suggested that different vibrotactile paradigms target different cortical 
mechanisms, and a number of recent studies have established links between somatosensory 
cortical function and measurable aspects of behavior. The relationship between cortical 
mechanisms and sensory function is particularly relevant with respect to developmental disorders 
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in which altered inhibitory processing has been postulated, such as in ASD and ADHD. In this 
study, a vibrotactile battery consisting of nine tasks (incorporating reaction time, detection 
threshold, and amplitude- and frequency discrimination) was applied to a cohort of healthy adults 
and a cohort of typically developing children to assess the feasibility of such a vibrotactile battery 
in both cohorts, and the performance between children and adults was compared. These results 
showed that children and adults were both able to perform these tasks with a similar performance, 
although the children were sUghtly less sensitive in frequency discrimination. Performance within 
different task-groups clustered together in adults, providing further evidence that these tasks tap 
into different cortical mechanisms, which is also discussed. This clustering was not observed in 
children, which may be potentially indicative of development and a greater variability. In 
conclusion, in this study, we showed that both children and adults were able to perform an 
extensive vibrotactile battery, and we showed the feasibility of applying this battery to other (e.g., 
neurodevelopmental) cohorts to probe different cortical mechanisms. 
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1. Introduction 

The cortical dynamics of somatosensory processing can be investigated using vibrotactile 
psychophysics. It has been suggested that different vibrotactile paradigms target different 
cortical mechanisms, and a number of recent studies have established links between 
somatosensory cortical function and measurable aspects of behavior (Hernandez et al., 2000; 
Puts et al., 201 1; Romo et al., 2003). However, links between GABAergic inhibitory 
neurotransmission and behavioral measures are less well understood. GABAergic inhibition 
is important in shaping the neuronal response to sensory stimulation (Alloway and Burton, 
1986; Dykes et al., 1984; Juliano et al., 1989), and most vibrotactile tasks rely in part on 
cortical GABAergic inhibitory mechanisms (Tommerdahl et al., 2010). Recent 
developments have made it possible to measure neurotransmitter concentration 
noninvasively in humans and correlate these concentrations with measures of tactile 
sensitivity (e.g.. Puts et al., 201 1). 

The relationship between GABA and sensory function is particularly relevant with respect to 
developmental disorders in which altered GABAergic processing has been postulated. For 
example, in Autism Spectrum Disorder (ASD), abnormal cortical structure (Casanova, 2004) 
and sensory processing (Blakemore et al., 2006; Tommerdahl et al., 2008b) have been linked 
to GABAergic processing, and GABA-system genetic variants have been proposed as 
models for ASD (e.g., (DeLorey, 2005)). In Tourette syndrome, both an altered density of 
GABAergic neurons (Kalanithi et al., 2005) and sensory impairments have been described 
(Belluscio et al., 2011; Miguel et al., 2000), and GABA gene markers correlate with tic 
severity (Tian et al., 201 1). Finally, GABA reductions have been shown in attention-deficit 
hyperactivity disorder (ADHD) (Edden et al., 2012), and impaired inhibition during cortical 
stimulation suggests reduced abnormal GABA interneuron activity (Gilbert et al., 201 1). 
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Thus, understanding the differences in sensory processing between groups may allow for a 
better understanding of cortical (dys)function in health and disease. 

In this study, we present a battery of vibrotactile tasks that targeted different aspects of 
cortical function. We demonstrate their feasibility in healthy adults (HA) and typically 
developing children (TDC), a prerequisite for future clinical studies, and present normative 
results. We present these data in the context of previous work in the field (Puts et al., 201 1; 
Tannan et al., 2007a, 2007b; Tommerdahl et al., 2008b; Zhang et al., 2011) (Lee et al., 2012; 
Nelson et al., 2012; Nguyen et al., 2013a, 2013b; Rai et al., 2012) to compare the 
performance of children and adults and to investigate patterns of performance. A priori, we 
would expect absolute levels of performance to differ between the HA and TDC but the 
relationships between related tasks to be preserved. 

1 .1 . Overview of task groups 

1 .1 .1 . Reaction time — A simple reaction time experiment ('press when you feel the 
stimulus') is a straightforward task for naive participantsthat allows them to become 
familiarized with the vibrotactile stimulation. Reaction time has been closely linked to white 
matter structure (Kerchner et al., 2012; Tamnes et al., 2012) and GABA concentration 
(Stagg et al., 201 1) in healthy subjects. In addition, reaction time has been shown to be 
altered in developmental disorders (Debes et al., 2011; Schuerholz et al., 1996; Xiao et al., 
2012). Reaction time probes both attentional and sensorimotor components. 

1.1.2. Detection threshold — The static detection threshold task is a well-known 
diagnostic tool. An abnormal detection threshold has been used as an indicator of brain 
dysfunction (Belluscio et al., 2011; Nudo et al., 2000; Staines et al., 2002) and is dependent 
on both white matter structure (Mountcastle et al., 1972) and GABAergic mechanisms 
(DeLorey et al., 201 1; Tavassoli et al., 2012). In a static vibrotactile detection threshold 
experiment, the weakest detectable stimulus is typically determined in either a yes/no or a 
two-alternative forced-choice (2AFC) manner. In contrast, a dynamic vibrotactile detection 
threshold experiment consists of a stimulus that is increased until perceived (see Zhang et 
al., 2011). It is thought that pre-detection sub-threshold stimulation mainly activates local 
feed-forward inhibitory mechanisms (Blankenburg et al., 2003; Favorov and Kursun, 201 1; 
Middleton et al., 2012; Swadlow, 2003), which thereby raises the detection threshold. 
Comparing dynamic and static threshold measures probes this feed-forward inhibition. 

1.1.3. Amplitude discrimination — Discriminating between two stimuli that are 
simultaneously applied to adjacent digits engages lateral inhibition to separate the response 
functions of the cortical areas representing each stimulus. A repetitive or 'adapting' stimulus 
has been shown to sharpen this response function (Whitsel et al., 1989, 2003), either by 
improving signal-to-noise or spatial resolution. Behaviorally, Hollins and Goble (1993) have 
shown that single-digit amplitude discrimination is improved by a 5 s adapting stimulus 
prior to each trial. In a similar fashion, Tannan et al. (2007b) have shown that in a healthy 
population, dual-site amplitude discrimination is improved when each trial is preceded by 
dual-site adaptation but is diminished when each trial is preceded by adaptation on only one 
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of the digits. Interestingly, this effect of adaptation is absent in ASD (Tommerdahl et al., 
2007). 

1 .1 .4. Frequency discrimination — Discriminating the frequency of two sequentially 
applied stimuli relies upon temporal processing. McLaughlin and Juliano (2005) showed 
that frequencies were, at least in part, encoded by the periodic synchronous firing of 
neuronal ensembles in the primary somatosensory cortex (SI) and that applying a GAB A 
antagonist destroys this periodicity. We have previously shown that individual differences in 
frequency discrimination performance were correlated with GABA concentration in the 
sensorimotor cortex, as measured by edited MRS (Puts et al., 201 1). In contrast, when 
frequencies are applied simultaneously to adjacent digits, temporal synchronization between 
the cortical areas, mediated by GABAergic lateral inhibition, would be expected to disrupt 
the temporal and periodic encoding of each stimulus, thereby impairing discrimination (e.g., 
Tommerdahl et al., 2008b). 

2. Materials and methods 

2.1. Participants 

Two cohorts were tested on a tactile battery consisting of nine tasks. Thirteen healthy adults 
(aged 30.5 + 4.9 years old; 3 female) and 22 typically developing children (aged from 8 to 
12 years old; 2 female) participated in this study. All of the participants were right-handed, 
which was confirmed using the Edinburgh Handedness Inventory (Oldfield, 1971) in the 
TDC cohort and by oral report in the healthy adult cohort. All of the TDC were recruited as 
controls for ongoing studies of ASD and ADHD. In TDC, the Wechsler Intelligence Scale 
for Children Third or Fourth Edition (WISC-III/IV) was used to assess intellectual ability. 
Children with full-scale IQ scores below 80 were excluded from participation (average IQ 
1 14.5 + 11.6). All of the children in the TDC cohort were free of criteria for psychiatric 
disorders as tested by the Diagnostic Interview for Children and Adolescents-Fourth Edition 
(DICA-IV), and none of the children in the TDC cohort were prescribed psychoactive 
medications. Informed consent was obtained from adult subjects and a parent of each child 
(who themselves assented to testing), under the approval of the Kennedy Krieger Institute 
and The Johns Hopkins School of Medicine Institutional Review Boards. 

2.2. Stimulus delivery 

A CM4 four-digit tactile stimulator (Cortical Metrics) was used for stimulation (Holden et 
al., 2012). All of the stimuli were delivered to the glabrous skin of the left digit 2 (LD2) and 
digit 3 (LD3) using a cylindrical probe (5 mm in diameter), and all stimuli were in the flutter 
range (25-50 Hz). Visual feedback, task responses, and data collection was performed using 
an Acer Onebook Netbook computer, running CM4 software (Holden et al., 2012). 

2.3. Experimental design 

The vibrotactile testing battery consisted of nine tasks, as shown in schematic form in Fig. 1. 
Prior to each task, the participants had to correctly respond to three consecutive practice 
trials to proceed, to confirm that the subject understood the instructions. Feedback was given 
during the practice trials but not during the task trials. In all tasks, stimulus delivery was 
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pseudo-randomized between LD2 and LD3. The response was obtained via a mouse-click 
using the participant's right hand. The left mouse button corresponded to LD3 and the right 
mouse button to LD2. All of the data were visually inspected prior to analysis. 

2.3.1. Reaction time: simple (sRT) and choice (cRT) reaction time — A 

suprathreshold stimulus (frequency 25 Hz, amplitude 300 [im, duration 40 ms) was delivered 
on LD2 or LD3, and the participants were asked to respond as quickly as possible when they 
felt the stimulus. In the sRT task, a mouse-click was sufficient, whereas in the cRT task, the 
participants additionally had to determine on which finger they felt the stimulus (inter-trial 
interval (ITI) 3 s; 20 trials). For each individual, the reaction times (for correct trials only in 
the cRT task) were sorted in ascending order, and the mean of the median 6 values was 
obtained as the mean RT. 

2.3.2. Detection thireshiold: static (sD) and dynamic (dD) detection thireshiold — 

In the sD task, a supra-threshold stimulus (frequency 25 Hz, starting amplitude 25 |jm, 
duration 500 ms) was delivered to either digit and the participants were asked to respond on 
which finger they felt the stimulus. A 1 up/1 down tracking paradigm (stimulus amplitude 
was decreased for a correct answer and increased for an incorrect answer) was used for the 
first 10 trials and a 2 up/1 down (two correct answers were necessary for a reduction in test 
amplitude) was used for the remainder of the task (ITI 5 s; 24 trials). The sD threshold was 
obtained as the mean amplitude of the final four trials, and the amplitude was determined for 
the twenty-fifth trial. In the dD task, after a variable delay (0-2500 ms), a 25 Hz stimulus 
increased from zero amplitude (rate of amplitude increase 2 |jm/s). The participants were 
asked to respond as soon as they felt the stimulus and to indicate the finger on which the 
stimulus was felt (ITI 10 s; 7 trials). The DD threshold was obtained as the mean stimulus 
amplitude at the time of pressing the button, across all correct trials. 

2.3.3. Amplitude discrimination threshold with no adaptation (nAD), dual-site 
adaptation (dAD) and single-site adaptation (sAD) — The amplitude discrimination 
tasks have been previously described (Tannan et al., 2007b; Tommerdahl et al., 2008a). In 
the nAD task, the participants were asked to choose which of the two simultaneously 
deUvered stimuli had the higher amplitude (25 Hz; 500 ms; Standard stimulus amplitude: 
100 [im; initial comparison stimulus amplitude: 200 |im). A 1 up/1 down tracking paradigm 
(comparison stimulus amplitude was decreased for a correct answer and increased for a 
wrong answer) was used for the first 10 trials and a 2 up/1 down (two correct answers were 
necessary for a reduction in comparison stimulus amplitude) was used for the remainder of 
the task (ITI 5 s; 20 trials). In the dAD condition, each trial was preceded by dual-site- 
delivered adapting stimuli (25 Hz; duration 1 s, amplitude 100 (im) and in the sAD task, 
each trial was preceded by a single-site-delivered adapting stimulus (duration 1 s, amplitude 
100 pm). Amplitude discrimination thresholds were obtained as the mean amplitude of the 
final four trials, and the amplitude was determined for the twenty-first trial. 

2.3.4. Frequency discrimination threshold: sequential (sqFD) and 
simultaneous (smFD) — In the sqFD task, stimuU (500 ms; 200 |jm) were deUvered to 
LD2 and LD3 sequentially (inter-stimulus interval 500 ms; pseudo-random location). In the 
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smFD task, the two stimuli were delivered simultaneously to both LD2 and LD3 
(pseudorandom location). One finger always received the standard stimulus (30 Hz) and the 
other the comparison stimulus (initial frequency 40 Hz). In both conditions, the participants 
were asked which finger received the higher frequency stimulus. The 1 up/1 down tracking 
paradigm (comparison stimulus frequency was decreased for a correct answer and increased 
for a wrong answer) was used for the first 10 trials and the 2 up/1 down (two correct 
answers were necessary for a reduction in comparison stimulus frequency) was used for the 
remainder of the task (ITI 5 s; 20 trials). Frequency discrimination thresholds were obtained 
as the mean of the frequency of the final four trials, and the frequency was determined for 
the twenty-first trial. Previous frequency discrimination studies have shown that the 
perceived intensity varies as a function of frequency as well as intensity (LaMotte and 
Mountcastle, 1975; Verrillo and Capraro, 1975). However, Harris et al. (2001) previously 
reported that the "subjects' accuracy at comparing frequency was not affected by shifts in 
vibration amplitude that causes the two vibrations to have equivalent intensity (i.e., by 
increasing the amplitude of a lower frequency)". In this study, the ampUtude was constant 
for both the standard and comparison stimuli, and the order of higher/lower was randomized 
across digits. 

2.4. Analysis 

The participants' data for an individual task were excluded when it was reported - orally by 
the experimenter- that the participant did not execute the task properly and showed poor 
behavioral compliance (e.g., pressing buttons as quickly as possible without regard for the 
stimulus and task), or when inspection of the tracking-profile showed large deviations in 
stimulus value over the last five trials (greater than four times the starting value, divided by 
the number of trials). Initial analysis focused on comparisons between related tasks (i.e., 
paired f-test between sRT and cRT; paired f-test between sD and dD; paired f-test between 
dAD and sAD; ANOVA of the three AD tasks; paired f-test of smFD and sqFD) keeping 
HA and TDC separate. The correlation matrices were calculated for all nine tasks for the HA 
and TDC groups independently. A clustering den-drogram was calculated from the 
normalized correlation matrix to investigate the relationships between the different tasks. 

3. Results 

One participant was fully excluded from the TDC cohort due to poor execution of all tasks, 
based on an oral report by the experimenter and confirmed visual inspection of data. 

3.1. Reaction time 

The mean RT for the sRT task and cRT tasks were 227.03 + 71.61 ms and 41 1.2 + 71.07 ms, 
respectively, for the HA group (an average increase in RT of 90% between the sRT and 
cRT) and 320.72 + 84.51 and 640.6803 + 191.60 ms for the TDC (an average increase in RT 
of 108% between the sRT and cRT), as shown in Fig. 2a. As expected, the mean RT was 
slower for the cRT task than for the sRT tasks for both groups (paired f-test p < 0.0001). The 
RT was significantly slower in both tasks for TDC compared to HA (p < 0.01 for both 
tasks), although the increase in RT was not. 
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3.2. Detection threshold 

For the HA group, the mean sD threshold was 4.84 +1.3 [im and the mean dD threshold was 
8.4 + 2.05 [im. For the TDC group, the mean sD threshold was 5.75 + 2.28 [im and the mean 
dD threshold was 8.68 ± 2.35 |im. As the stimulus amplitude continued to increase between 
perception and response, the reaction time component in the dD task increased its threshold 
values; thus, the dD threshold for each individual was corrected using their mean cRT and 
the rate of amplitude increase. This resulted in corrected dD thresholds of 7.44 + 2.16 pm for 
the HA group and 7.42 + 2.32 pm for the TDC group. The sD and corrected dD are shown in 
Fig. 2b. The dD threshold was significantly greater than the sD threshold in both groups (p < 
0.001). The results from one additional participant in the TDC group were excluded due to 
poor execution of the task. There were no differences in task performance or difference in 
sD-dD difference between the two cohorts (p > 0.5). 

3.3. Amplitude discrimination 

In the HA group, the mean AD threshold was 46.15 + 18.5 |jm without adaptation, 34.38 + 
18.57 |jm with dual-site adaptation (an average increase in performance of 21%) and 58.08 + 
21.84 |jm with single-site adaptation (an average decrease in performance of 36%), as shown 
in Fig. 2c. A one-way ANOVA showed a significant difference between the task 
performance in the three adaptation conditions, and post hoc paired f-tests showed that the 
dAD and sAD thresholds were significantly different (p = 0.013) and that the differences 
between dAD and nAD and between sAD and nAD were close to threshold (p = 0.055, and 
p = 0.0494, respectively). 

For the TDC group, the mean AD threshold was 49.21 + 29.98 pm without adaptation, 39.41 
+ 22.20 |im with dual-site adaptation (an average increase in performance of 20%) and 65.81 
+ 36.22 |jm with single-site adaptation (an average decrease in performance of 34%). A one- 
way ANOVA showed a significant difference (p < 0.05) between the conditions, and post 
hoc paired f-tests showed that both nAD and dAD thresholds were significantly different 
from the sAD thresholds (p = 0.010, and p = 0.001, respectively) but that they did not differ 
from one another (p = 0.7). Four TDC participants were excluded from the amplitude 
discrimination task (one for nAD, three for dAD) due to poor execution of the task and 
improper tracking. Moreover, there were no significant differences in AD performance 
between the two cohorts (p = 0.08 for nAD, p > 0.5 for dAD and sAD). 

3.4. Frequency discrimination 

For the HA group, the mean frequency discrimination threshold was 5.4 + 2.4 Hz in the 
sequential condition and 10.2 + 3.5 Hz in the simultaneous condition (significantly different 
atp < 0.0001), as shown in Fig. 2d. For the TDC group, the mean frequency discrimination 
threshold was 7.68 + 2.31 Hz in the sequential condition and 9.31 + 2.61 Hz in the 
simultaneous condition (not significantly different alp < 0.05). In the TDC group, one 
participant was excluded for both tasks due to poor execution, three were excluded for the 
sequential condition due to poor tracking, and one was excluded from the simultaneous task 
due to poor tracking. The TDC performed significantly worse than the adults in the 
sequential, but not in the simultaneous FD task (p < 0.05, p > 0.4, respectively). 
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3.5. Correlation analysis of all tasks 

The correlation matrix for the HA group, shown in Fig. 3 a, demonstrates particularly strong 
relationships between the two RT tasks, between the FD tasks, and between the three 
amplitude discrimination tasks. Furthermore, the amplitude discrimination following single- 
site adaptation was negatively correlated with corrected dynamic detection threshold (R = 
-0.83). There was a correlation between the dD threshold and cRT, which decreased after 
correcting for reaction time, suggesting an important role of reaction time in the dD 
threshold {R = 0.5 and 0.38, respectively). The clustering dendrogram in the analysis 
depicted in Fig. 3b showed that the related tasks clustered together (RT, R = 0.78; AD, R = 
0.35-0.46; FD, R = 0.42), although the dynamic detection threshold (both corrected and 
uncorrected dD) clustered to some extent with the RT tasks. The correlation matrix for the 
TDC group, as shown in Fig. 3a, appeared to show more, but weaker, correlations. Similar 
to the HA group, the RT tasks were correlated with each other as well (R = 0.49) and with 
both detection threshold tasks (R = 0.38). Consistent with the HA group, the dynamic 
detection threshold was negatively correlated with amplitude discrimination following 
single-site adaptation {R = -0.51). However, the correlations among the amplitude 
discrimination tasks and among the frequency discrimination tasks as shown in the HA 
group were absent in the TDC group, and the dendrogram appeared uninformative. 

4. Discussion 

We have presented vibrotactile behavioral data on a battery of tasks that can be collected 
from cohorts of healthy adults (HA) and typically developing children (TDC) in a total 
testing time of approximately 30 min. Both HA and TDC of 8-12 years of age were able to 
perform these tasks. Moreover, we have shown not only that TDC are able to perform these 
tasks but also that they show the same patterns of performance as healthy adults (the effect 
of feed-forward inhibition on detection threshold, the effect of adaptation on amplitude 
discrimination and the effect of synchronous stimulation on frequency discrimination), 
which suggested that the mechanisms underlying these tasks were similar between healthy 
adults and healthy children. It should be noted that three adult female participants performed 
the tasks. It is possible that menstrual cycle has an effect on tactile sensitivity, although the 
evidence is inconsistent (Bajaj et al., 2001, 2002; Drobek et al., 2002). All of the participants 
performed the tasks within a single session, and any differences in sensitivity due to 
hormonal effects were expected to be reflected within the variance of the data. 

In both HA and TDC, as expected, the choice reaction times were longer than the simple 
reaction times, and the two results were highly correlated across individuals within both 
groups. Children performed significantly more slowly than adults in both reaction time 
tasks, consistent with previous studies showing a U-shaped relationship between age and 
reaction time (Williams et al., 2005). 

The static detection threshold was slightly higher for TDC, which was consistent with 
previous findings (Bernstein et al., 1986; Guclu and Oztek, 2007), although it was not 
significant. There was no significant difference between groups in the dynamic detection 
threshold and no difference at all after correction for reaction time. It has previously been 
shown that the dD threshold was greater than sD threshold in HA and that adults become 
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worse at both tasks with age (Zhang et al., 201 1). To the best of our knowledge, this task has 
not previously been examined in children. The difference between the dynamic and static 
detection thresholds has been suggested to be related to feed-forward GABAergic inhibition. 
Favorov and Kursun (201 1) have suggested that layer IV is highly involved in computations 
of the feed-forward inhibitory drive, which is affected by prior stimulus information. 
Blankenburg et al. (2003) showed that a preliminary sub-threshold stimulus 30 ms prior to 
the detection trial increases the detection thresholds, and the authors proposed that this effect 
was due to cortical feed-forward inhibitory mechanisms, as inhibitory internem-ons have a 
lower spiking threshold than excitatory neurons and are therefore more strongly activated by 
sub threshold stimuli (Gil and Amitai, 1996). Blankenburg et al. (2003) discussed the 
possibility that this feed-forward inhibition might be protective against spurious activity in 
the cortex by decreasing the net cortical activity. Zhang et al. (201 1) suggested that sub- 
threshold stimulation in the dynamic threshold task drives inhibitory mechanisms, which 
was supported by the observation that the difference diminishes with age (Zhang et al., 
2011) as GABAergic inhibition declines. It also appears that this feed-forward processing 
was fully developed in the cohort examined (8-12 years old). 

In amplitude discrimination, dual-site adaptation tends to improve performance, whereas 
single-site adaptation worsens it, as expected from previous studies (Tannan et al., 2007 a,b; 
Zhang et al., 201 1). Our results were consistent with (although extrapolatory to) previous 
reports that amplitude discrimination performance and the effect of adaptation did not 
change with age (Zhang et al., 201 1). 

In frequency discrimination, both cohorts performed worse in the simultaneous frequency 
discrimination task. In fact, most subjects had significant difficulty in advancing beyond the 
practice trials and reported being unable to perform the task (some children were unable to 
perform either task, supported by an oral report from both children and experimenters, 
potentially due to frequency as an abstract concept). The mean simultaneous FD threshold 
was not significantly different from the initial difference of 10 Hz. It appeared that whereas 
sequentially applied stimuli can be distinguished on the basis of frequency, simultaneously 
applied stimuli cannot, or at least not better than 10 Hz. It is possible that synchronization 
between cortical representations of the two stimuli, mediated by lateral inhibitory pathways, 
reduced the perceptual separation of signals and impaired task performance. The detrimental 
effect of cortical synchronization between digits on stimulus separation has been previously 
described (Tommerdahl et al., 2008), with reference to temporal order judgment. In that 
study, negative effects of synchronization were not observed in participants with autism, and 
the authors suggested this phenomenon could be due to reduced GABAergic local 
connectivity. 

The percepts of frequency and amplitude are not independent, and higher frequencies tend to 
be perceived as having higher amplitude (LaMotte and Mountcastle, 1975; Verrillo and 
Capraro, 1975). Some studies (Goble and HoUins, 1994) have used frequency-amplitude 
matching to remove amplitude as a potential driver of frequency discrimination performance 
(also; Goble and Hollins, 1994). The aim of this study was to develop a short battery of 
tasks. Thus, we used a frequency discrimination task that involved physically equal rather 
than perceptually matched stimulus amplitudes, as has been done by Harris and colleagues 
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(e.g., (Harris et al., 2001, 2006)). Harris et al. (2001) found that amplitude matching did not 
reduce frequency discrimination performance, and interestingly, our current results showed 
no significant correlation between amplitude and frequency discrimination performance. 
One potential improvement on physically matched amplitudes would be to introduce an 
amplitude jitter. 

Interestingly, the correlational and clustering analysis in the HA group tended to sort tasks 
by property. Thus, the subjects who performed well at the simple reaction time task tended 
to perform well at the choice reaction time; subjects who performed well at the static 
detection threshold task also tended to perform well at the dynamic task; etc. This outcome 
was perhaps not surprising; however, the fact that the two frequency discrimination tasks 
clustered suggested that it was not true to simply imply that the subjects 'cannot perform the 
task.' The modulus mean between-task correlation coefficient was 0.23 for HA and 0.26 for 
children (compared to the maximum correlation value of 0.78 in HA), suggesting that it was 
not true that these tasks were equivalent and that the subjects did not perform 'well' or 
'badly' in equal measure across tasks. This finding was the main value of performing a 
battery of tasks such as the one presented. In general, the clustering of similar tasks 
suggested that measuring performance with respect to different properties of vibrotactile 
stimuli targeted different aspects of cortical processing. 

The correlational analysis did reveal some links between task groupings. The relationship 
between detection threshold and reaction time was not expected, although both tasks were 
related to white matter integrity (Tamnes et al., 2012; Kerchner et al., 2012). The negative 
correlation between detection threshold and amplitude discrimination (non-adapted and with 
single-site adaptation) makes sense in terms of cortical inhibition; subjects with greater 
levels of inhibition will have higher detection thresholds but would be expected to show 
better discrimination. Interestingly, this relationship was not maintained for dual-site 
adaptation, although this result could be due to a bottom effect and there might not be room 
for improvement in the dual-site adaptation condition. In addition, in neither HA nor TDC 
were frequency and amplitude discrimination correlated, which may provide additional 
evidence that ampUtude information is not used in the frequency discrimination task. 

Correlation and dendrogram analyses revealed differences in the way the domains are 
related between HA and TDC. However, neither the frequency tasks nor the amplitude 
discrimination tasks appeared to be correlated in the TDC, while they were in HA. 
Differences in task relationships between HA and TDC may be indicative of development 
(neither cutaneous nerves nor spinal cord mature fully until after puberty (e.g., (Allison et 
al., 1984; Sato et al., 1977)), sensorimotor development, or the development of attentional 
control. It is possible that some children have more difficulty with tasks than others, thereby 
increasing variability and masking correlational relationships. 

4.1. Limitations 

While it is beneficial to present a battery of tasks to investigate a number of different 
processes, these methods do have some limitations. The aim of this battery was that it could 
be performed within 30 min, which makes it suitable for a naive cohort as well as for 
pediatric populations. However, most behavioral tasks described in the psychophysics 
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literature are lengthy (in terms of trial numbers and task repeats), and it remains unclear 
what effect our shorter protocols have on the accuracy of the measurements. However, the 
brevity of testing does allow for a much greater number of participants to be tested at the 
same time, partially offsetting any loss of statistical power. 

The task duration was reduced by starting the tasks with relatively difficult initial settings, 
simultaneous presentation of stimuli to two digits (i.e., smFD is shorter than sqFD), and the 
randomization of the order of stimuli and parameters not being tested (e.g., amplitude was 
pseudo-random in the frequency discrimination tasks to reduce discrimination on the basis 
of amplitude, without having to perform frequency-amplitude correction). Visual inspection 
of the tracking curves showed that the majority of the participants reached a plateau. The 
threshold was obtained as the mean of the final five trials, and the average coefficient of 
variation in adults for all tasks was less than 10% of the average threshold value, indicating 
a small variability of the last five trials. In addition, the results shown by our naive cohort 
compared well with previous findings on these tasks. Testing a pediatric cohort can be 
challenging, and increasing the number of trials, while potentially increasing the SNR of the 
measurement, would be expected to adversely affect compliance and threshold 
measurements. Testing the reproducibility of these measurements within naive cohorts is 
problematic because a number of different studies have shown effects of perceptual learning 
in these tasks. However, the strong within-task-group correlations observed are at least 
circumstantial evidence for good within-task reproducibility. 

In conclusion, we have presented a 30 min tactile behavioral battery that probes a number of 
different cortical mechanisms and that is easily applied to adults and children as young as 8 
years old. 

Acknowledgments 

NAJP is funded by an Autism Speaks translational postdoctoral fellowsliip. ELW is supported by the Organization 
for Autism Research. This work was further supported by NIH; P41 EBO 15909 and NIH/NINDS: 2 ROl 
NS048527-08; 2 ROl MH078160 and The Johns Hopkins University School of Medicine Institute for Clinical and 
Translational Research National Institutes of Health/National Center for Research Resources Clinical and 
Translational Science Award program ULl RR025005. We would like to thank Taylor Koriakin, Carrie Nettles and 
Rebecca Buhlman for their help in acquiiing the TDC data. 

Abbreviations 



HA 


healthy adults 


TDC 


typically developing children 


LD2/LD3 


left digit 2 and left digit 3 


sRT 


simple reaction time task 


CRT 


choice reaction time task 


sD 


static detection threshold task 


dD 


dynamic detection threshold task 


nAD 


amplitude discrimination - no adaptation 
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dAD 



amplitude discrimination - dual-site adaptation 



sAD 



amplitude discrimination - single-site adaptation 



sqFD 



sequential frequency discrimination 



smFD 



simultaneous frequency discrimination 



ISI 



interstimulus interval 



m 



intertrial interval 



References 



Allison T, Hume AL, Wood CC, Goff WR. Developmental and aging changes in somatosensory, 
auditory and visual evoked potentials. Electroencephalogr Clin Neurophysiol. 1984; 58:14-24. 
[PubMed: 6203699] 

AUoway KD, Burton H. BicucuUine-induced alterations in neuronal responses to controlled tactile 

stimuli in the second somatosensory cortex of the cat: a microiontophoretic study. Somatosens Res. 

1986; 3:197-211. [PubMed: 2875510] 
Bajaj P, Arendt-Nielsen L, Bajaj P, Madsen H. Sensory changes during the ovulatory phase of the 

menstrual cycle in healthy women. Eur J Pain. 2001 ; 5: 135-44. [PubMed: 1 1465979] 
Bajaj P, Bajaj P, Madsen H, Arendt-Nielsen LA. Comparison of modality-specific somatosensory 

changes during menstruation in dysmenorrheic and nondysmenorrheic women. Clin J Pain. 2002; 

18:180-90. [PubMed: 12048420] 
Belluscio BA, Jin L, Walters V, Lee TH, Hallett M. Sensory sensitivity to external stimuli in Tourette 

syndrome patients. Mov Disord. 2011; 26:2538-43. [PubMed: 22038938] 
Bernstein LE, Schechter MB, Goldstein MH Jr. Child and adult vibrotactile thresholds for sinusoidal 

and pulsatile stimuh. J Acoust Soc Am. 1986; 80: 118-23. [PubMed: 3745657] 
Blakemore SJ, Tavassoli T, Calo S, Thomas RM, Catmur C, Frith U, Haggard P. Tactile sensitivity in 

Asperger syndrome. Brain Cogn. 2006; 61:5-13. [PubMed: 16500009] 
Blankenburg F, Taskin B, Ruben J, Moosmann M, Ritter P, Curio G, Villringer A. Imperceptible 

stimuli and sensory processing impediment. Science. 2003; 299:1864. [PubMed: 12649475] 
Casanova MF. White matter volume increase and minicolumns in autism. Ann Neurol. 2004; 56:453. 

[author reply 4]. [PubMed: 15349878] 
Debes NM, Hansen A, Skov L, Larsson HA. functional magnetic resonance imaging study of a large 
clinical cohort of children with Tourette syndrome. J Child Neurol. 201 1; 26:560-9. [PubMed: 
21464239] 

DeLorey TM. GABRB3 gene deficient mice: a potential model of autism spectrum disorder. Int Rev 

Neurobiol. 2005; 71:359-82. [PubMed: 16512358] 
DeLorey TM, Sahbaie P, Hashemi E, Li WW, Salehi A, Clark DJ. Somatosensory and sensorimotor 

consequences associated with the heterozygous disruption of the autism candidate gene, Gabrb3. 

Behav Brain Res. 2011; 216:36-45. [PubMed: 20699105] 
Drobek W, Schoenaers J, De Laat A. Hormone-dependent fluctuations of pressure pain threshold and 

tactile threshold of the temporalis and masseter muscle. J Oral Rehabil. 2002; 29:1042-51. 

[PubMed: 12453257] 

Dykes RW, Landry P, Metherate R, Hicks TP. Functional role of GABA in cat primary somatosensory 
cortex: shaping receptive fields of cortical neurons. J Neurophysiol. 1984; 52:1066-93. [PubMed: 
6151590] 

Edden RA, Crocetti D, Zhu H, Gilbert DL, Mostofsky SH. Reduced GABA concentration in attention- 
deficit/hyperactivity disorder. Arch Gen Psychiatry. 2012; 69:750-3. [PubMed: 22752239] 

Favorov OV, Kursun O. Neocortical layer 4 as a pluripotent function linearizer. J Neurophysiol. 201 1; 
105:1342-60. [PubMed: 21248059] 



J Neurosci Methods. Author manuscript; available in PMC 2014 August 15. 



Puts et al. 



Page 13 



Gil Z, Amitai Y. Properties of convergent thalamocortical and intracortical synaptic potentials in 
single neurons of neocortex. J Neurosci. 1996; 16:6567-78. [PubMed: 8815933] 

Gilbert DL, Isaacs KM, Augusta M, Macneil LK, Mostofsky SH. Motor cortex inhibition: a marker of 
ADHD behavior and motor development in children. Neurology. 201 1; 76:615-21. [PubMed: 
21321335] 

Goble AK, HoUins M. Vibrotactile adaptation enhances frequency discrimination. J Acoust Soc Am. 

1994; 96:771-80. [PubMed: 7930078] 
Guclu B, Oztek C. Tactile sensitivity of children: effects of frequency, masking, and the non-Pacinian 

I psychophysical channel. J Exp Child Psychol. 2007; 98:113-30. [PubMed: 17597140] 
Harris JA, Arabzadeh E, Fairhall AL, Benito C, Diamond ME. Factors affecting frequency 

discrimination of vibrotactile stimuli: implications for cortical encoding. PLoS ONE. 2006; 

l:el00. [PubMed: 17183633] 
Harris JA, Harris IM, Diamond ME. The topography of tactile learning in humans. J Neurosci. 2001; 

21:1056-61. [PubMed: 11157091] 
Hernandez A, Zainos A, Romo R. Neuronal correlates of sensory discrimination in the somatosensory 

cortex. Proc Natl Acad Sci USA. 2000; 97:6191-6. [PubMed: 1081 1922] 
Holden JK, Nguyen RH, Francisco EM, Zhang Z, Dennis RG, Tommerdahl M. A novel device for the 

study of somatosensory information processing. J Neurosci Methods. 2012; 204:215-20. 

[PubMed: 22155443] 

Hollins M, Goble AK. Vibrotactile adaptation enhances amplitude discrimination. J Acoust Soc Am. 

1993; 93:418-24. [PubMed: 8423258] 
Juliano SL, Whitsel BL, Tommerdahl M, Cheema SS. Determinants of patchy metabolic labeling in 

the somatosensory cortex of cats: a possible role for intrinsic inhibitory circuitry. J Neurosci. 1989; 

9:1-12. [PubMed: 2913199] 
Kalanithi PS, Zheng W, Kataoka Y, DiFiglia M, Grantz H, Saper CB, Schwartz ML, Leckman JF, 

Vaccarino FM. Altered parvalbumin-positive neuron distribution in basal ganglia of individuals 

with Tourette syndrome. Proc Natl Acad Sci USA. 2005; 102:13307-12. [PubMed: 16131542] 
Kerchner GA, Racine CA, Hale S, Wilheim R, Laluz V, Miller BL, Kramer JH. Cognitive processing 

speed in older adults: relationship with white matter integrity. PLoS ONE. 2012; 7:e50425. 

[PubMed: 23185621] 

LaMotte RH, Mountcastle VB. Capacities of humans and monkeys to discriminate vibratory stimuli of 

different frequency and amplitude: a correlation between neural events and psychological 

measurements. J Neurophysiol. 1975;38:539-59. [PubMed: 1127456] 
Lee J, Woo J, Favorov OV, Tommerdahl M, Lee CJ, Whitsel BL. Columnar distribution of activity 

dependent gabaergic depolarization in sensorimotor cortical neurons. Mol Brain. 2012; 5:33. 

[PubMed: 23006518] 

McLaughlin DF, Juliano SL. Disruption of layer 4 development alters laminar processing in fen'et 

somatosensory cortex. Cereb Cortex. 2005; 15:1791-803. [PubMed: 15772374] 
Middleton JW, Omar C, Doiron B, Simons DJ. Neural correlation is stimulus modulated by feed 

forward inhibitory circuitry. J Neurosci. 2012; 32:506-18. [PubMed: 22238086] 
Miguel EC, do Rosario-Campos MC, Prado HS, do Valle R, Ranch SL, Coffey BJ, Baer L, Savage 

CR, O' Sullivan RL, Jenike MA, Leckman JF. Sensory phenomena in obsessive-compulsive 

disorder and Tourette's disorder. J Clin Psychiatry. 2000; 61:150-6. quiz 7. [PubMed: 10732667] 
Mountcastle VB, LaMotte RH, Carli G. Detection thresholds for stimuli in humans and monkeys: 

comparison with threshold events in mechanoreceptive afferent nerve fibers innervating the 

monkey hand. J Neurophysiol. 1972; 35:122-36. [PubMed: 4621505] 
Nelson AJ, Premji A, Rai N, Hoque T, Tommerdahl M, Chen R. Dopamine alters tactile perception in 

Parkinson's disease. Can J Neurol Sci. 2012; 39:52-7. [PubMed: 22384496] 
Nguyen RH, Ford S, Calhoun AH, Holden JK, Gracely RH, Tommerdahl M. Neurosensory 

assessments of migraine. Brain Res. 2013a 
Nguyen RH, Gillen C, Garbutt JC, Kampov-Polevoi A, Holden JK, Francisco EM, Tommerdahl M. 

Centrally-mediated sensory information processing is impacted with increased alcohol 

consumption in college-aged individuals. Brain Res. 2013b; 1492:53-62. [PubMed: 23178333] 



J Neurosci Methods. Author manuscript; available in PMC 2014 August 15. 



Puts et al. 



Page 14 



Nudo RJ, Friel KM, Delia SW. Role of sensory deficits in motor impairments after injury to primary 

motor cortex. Neuropharmacology. 2000; 39:733-42. [PubMed: 10699440] 
Oldfield RC. The assessment and analysis of handedness: the Edinburgh inventory. Neuropsychologia. 

1971; 9:97-113. [PubMed: 5146491] 
Puts NAJ, Edden RAE, Evans CJ, McGlone F, McGonigle DJ. Regionally specific human GABA 

concentration correlates with tactile discrimination thresholds. J Neurosci. 2011; 31:16556-60. 

[PubMed: 22090482] 

Rai N, Premji A, Tommerdahl M, Nelson AJ. Continuous theta-burst rTMS over primary 

somatosensory cortex modulates tactile perception on the hand. Clin Neurophysiol. 2012; 

123:1226-33. [PubMed: 22100859] 
Romo R, Hernandez A, Zainos A, Salinas E. Correlated neuronal discharges that increase coding 

efficiency during perceptual discrimination. Neuron. 2003; 38:649-57. [PubMed: 12765615] 
Sato S, Ogihara Y, Hiraga K, Nishijima A, Hidano A. Fine structure of unmyelinated nerves in 

neonatal skin. J Cutan Pathol. 1977; 4:1-8. [PubMed: 886028] 
Schuerholz LJ, Baumgardner TL, Singer HS, Reiss AL, Denckla MB. Neuropsycholo-gical status of 

children with Tourette's syndrome with and without attention deficit hyperactivity disorder. 

Neurology. 1996; 46:958-65. [PubMed: 8780072] 
Stagg CJ, Bestmann S, Constantinescu AO, Moreno LM, AUman C, Mekle R, Woolrich M, Near J, 

Johansen-Berg H, Rothwell JC. Relationship between physiological measures of excitability and 

levels of glutamate and GABA in the human motor cortex. J Physiol. 201 1; 589:5845-55. 

[PubMed: 22005678] 

Staines WR, Black SE, Graham SJ, Mcllroy WE. Somatosensory gating and recovery from stroke 

involving the thalamus. Stroke. 2002; 33:2642-51. [PubMed: 1241 1655] 
Swadlow HA. Fast-spike interneurons and feedforward inhibition in awake sensory neocortex. Cereb 

Cortex. 2003; 13:25-32. [PubMed: 12466212] 
Tamnes CK, Fjell AM, Westlye LT, Ostby Y, Walhovd KB. Becoming consistent: developmental 

reductions in intraindividual variability in reaction time are related to white matter integrity. J 

Neurosci. 2012; 32:972-82. [PubMed: 22262895] 
Tannan V, Holden JK, Tommerdahl M, Favorov OV, Zachek M. Effects of stimulus-driven 

synchronization on sensory perception. Behav Brain Funct. 2007a; 3:61. [PubMed: 18053216] 
Tannan V, Simons S, Dennis RG, Tommerdahl M. Effects of adaptation on the capacity to 

differentiate simultaneously delivered dual-site vibrotactile stimuli. Brain Res. 2007b; 1 186: 164- 

70. [PubMed: 18005946] 
Tavassoli T, Auyeung B, Murphy LC, Baron-Cohen S, Chakrabarti B. Variation in the autism 

candidate gene GABRB3 modulates tactile sensitivity in typically developing children. Mol 

Autism. 2012; 3:6. [PubMed: 22769427] 
Tian Y, Gunther JR, Liao IH, Liu D, Ander BP, Stamova BS, Lit L, Jickling GC, Xu H, Zhan X, Sharp 

FR. GABA- and acetylcholine-related gene expression in blood correlate with tic severity and 

microarray evidence for alternative splicing in Tourette syndrome: a pilot study. Brain Res. 201 1; 

1381:228-36. [PubMed: 21241679] 
Tommerdahl M, Favorov OV, Whitsel BL. Dynamic representations of the somatosensory cortex. 

Neurosci Biobehav Rev. 2010; 34:160-70. [PubMed: 19732790] 
Tommerdahl M, Tannan V, Cascio CJ, Baranek GT, Whitsel BL. Vibrotactile adaptation fails to 

enhance spatial localization in adults with autism. Brain Res. 2007; 1154:116-23. [PubMed: 

17498672] 

Tommerdahl M, Tannan V, Francisco E, Zhang Z, Holden J. Vibrotactile amplitude discrimination 

capacity parallels magnitude changes in somatosensory cortex and follows Weber's Law. Exp 

Brain Res. 2008a; 191:49-56. [PubMed: 18651137] 
Tommerdahl M, Tannan V, Holden JK, Baranek GT. Absence of stimulus-driven synchronization 

effects on sensory perception in autism: Evidence for local underconnectivity? Behav Brain Funct. 

2008b; 4:19. [PubMed: 18435849] 
Verrillo RT, Capraro AJ. Effect of stimulus frequency on subjective vibrotactile magnitude functions. 

Perception. 1975; 17:91-6. 



J Neurosci Methods. Author manuscript; available in PMC 2014 August 15. 



Puts et al. 



Page 15 



Whitsel, BL.; Favorov, O.; Tommerdahl, M.; Diamond, M.; Juliano, SJ.; Kelly, D. Sensory processing 
in the mammalian brain: neural substrates and experimental strategies. Oxford University Press; 
1989. Dynamic processes govern the somatosensory cortical response to natural stimulation; p. 
84-116. 

Whitsel BL, Kelly EF, Quibrera M, Tommerdahl M, Li Y, Favorov OV, Xu M, Metz CB. Time- 
dependence of SI RA neuron response to cutaneous flutter stimulation. Somatosens Mot Res. 
2003; 20:45-69. [PubMed: 12745444] 

Williams BR, Hultsch DF, Strauss EH, Hunter MA, Tannock R. Inconsistency in reaction time across 
the life span. Neuropsychology. 2005; 19:88-96. [PubMed: 15656766] 

Xiao T, Xiao Z, Ke X, Hong S, Yang H, Su Y, Chu K, Xiao X, Shen J, Liu Y. Response inhibition 
impairment in high functioning autism and attention deficit hyper-activity disorder: evidence from 
near-infrared spectroscopy data. PLoS ONE. 2012; 7:e46569. [PubMed: 23056348] 

Zhang Z, Francisco EM, Holden JK, Dennis RG, Tommerdahl M. Somatosensory information 
processing in the aging population. Front Aging Neurosci. 2011; 3:18. [PubMed: 22163221] 



J Neurosci Methods. Author manuscript; available in PMC 2014 August 15. 



Puts et al. 



Page 16 



Reaction Time 



Simple (SRT) 

Digit 2 



-V- 



stimulus Button press 



Choice (CRT) 



Mr 



Digit 3 

b 

static (sD) 

Digit 2 



Empty 



n 



Stimuius Response 



Empty 



40 ms 



40 ms 



Detection Tlireshold 

Dynamic (dD) 



Respond when perceived 



Stimulus 



Response 




Digit 3 



Empty 



Empty 



No adaptation (nAD) 



Digit 2 



Amplitude Discrimination 

SIngle/Dual-site adaptation (dAD and sAD) 



Digit 3 



standard 
- ^ Response 




var comparison 



"Adaptation" 



''i !' ,", ;\ A ;■, ,"i I' 
I 1 V V V II V '.' V w 



Test-Stimulus 



standard 
■+* »^ Response 




var comparison 



Frequency Discrimination 

Simultaneous (smFD) Sequential (sqFD) 



Digit 2 



Digits 




.standard _ ^ Response 




comparison 




standard 



empty 



empty 




Response 



comparison 



Fig. 1. 

Vibrotactile testing battery, trial examples, (a) Simple (sRT) and choice (cRT) reaction time, 
(b) Static (sD) and dynamic (dD) detection threshold, (c) Amplitude discrimination without 
adaptation (nAD), with dual-site adaptation (dAD) and single-site adaptation (sAD). (d) 
Sequential (sqFD) and simultaneous (smFD) frequency discrimination. 
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Fig. 2. 

Individual results for all tasks, (a) Reaction time. RT was faster in the sRT task than in the 
cRT task in both HA and TDC {p < 0.001). The TDC were significandy slower (p < 0.01) 
than the HA. (b) Detection threshold. The sD was significantly lower than the dD in both 
HA and TDC {p < 0.001). There was no significant difference in the detection threshold 
between HA and TDC. (c) Amplitude discrimination. In HA, the sAD threshold was 
significantly worse than the dAD {p < 0.02) and close to significance from nAD {p = 0.0494, 
uncorrected for multiple comparisons). The dAD was close to being significantly different 
from nAD {p = 0.055). In TDC, the sAD was also significantly worse than the dAD and 
nAD ip < 0.02), but the dAD and nAD did not differ significantly. There were no 
differences between the cohorts, (d) Frequency discrimination. sqFD was significantly better 
than smFD in HA {p < 0.05), but not in TDC. *p < 0.05. Box plot whisker are 5th and 95th 
percentile, center of the box is the mean. 



J Neurosci Methods. Author manuscript; available in PMC 2014 August 15. 



Puts et al. 



Page 18 



Healthy Adults 




Typically Developing Children 




Fig. 3. 

Correlation matrix and cluster analysis, (a) In HA, three different task groupings (RT + DT; 
AD; FD) correlated with each other. Furthermore, the sAD was negatively correlated with 
the corrected dD threshold (R = -0.83). (b) Cluster analysis clustered different tasks-groups 
within separate branches, although dD clustered with RT to some extent, (c) In the TDC, the 
correlation matrix showed more, but weaker correlations. The RT tasks are correlated with 
each other (R = 0.49) and with both detection threshold tasks {R = 0.38). Consistent with the 
HA group, the dD threshold was negatively correlated with the sAD {R = -0.51). However, 
the correlations among the amplitude discrimination tasks and among the frequency 
discrimination tasks, as shown in the HA group, were absent in the TDC group, (d) No 
clustering could be observed in the TDC. 
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